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SUMMART 


The theoretical critical stress is derived for a simply supported 
curved rectangular panel in axial compression having a central chord- 
vise stiffener offering no torsional restraint. The results are 
presented in the form of computed curves and a table. 

Because a panel of moderate or large curvature ‘buckles in com- 
pression at a stress considerably belov the theoretical value, a 
method Is suggested to aid in determining the critical stress for use 
in design. 


INTRODUCTION 


A simplified method vas recently developed for determining the 
theoretical buckling stresses of unstiffened cylindrical shells undVr 
various loading conditions (references 1 and 2). In the present paper 
the theory is extended to include stiffened shells (appendix- B), and 
the particular case of curved rectangular panels In axial compression, 
reinforced by a centrally located chordvise stiffener of zero torsional 
stiffness. Is treated in detail (appendix C). Numerical results for 
this case are given In figures and a table. Because teste shov that 
unstiffened curved panels in compression buckle at a load considerably 
belov the theoretical (see reference 3)* & procedure is suggested to 
permit the estimation of the actual critical stress of a rectangular 
panel -with a central chordvise stiffener. 


RESULTS AND DISCUSSION 


Theoretical critical stress .— The effect of a central chordvise 
stiffener on the critical axial stress of a curved rectangular panel 
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is to increase the buckling stress. The increase is very small or 
zero for 3 % 0.7 (P is the ratio of the circumferential dimension 
to the axial dimension) or when Z (curvature parameter, see appendix A) 
is greater than the values given in the following table: 


3 

Z 

0.83 

20.0 

1.0 

14.0 

1.2 

11.2 

1.5 

9.3 

2.0 

7.6 

3.0 

6 . 3 


Only when P > 0.7 and Z is less than the corresponding value 
given in this table, is the buckling stress for the panel with a 
central chordwise stiffener appreciably greater than the buckling 
stress for the panel without a stiffener. The percentage increase 
that can be expected is shown for a number of cases in table 1. 

The critical stress for a curved plate with a central chordwise 
stiffener is conveniently computed by the use of the standard buckling 
equation for a curved plate (see equation in fig. l) with the buckling 
coefficient k x for the unstiffened plate increased by an amount 
due to the presence of the stiffener. The critical-stress coefficient' k x 
for an unstiffened curved plate of various aspect ratios may be obtained 
frcm figure 1. Figure l(a; gives the variation of k^. with aspect 
ratio for a flat plate (Z = 0)j whereas figure l(b) gives k x for 
curved plates. The increase in critical-axial— stress coefficient Ak^ 
due to the presence of the stiffener is given in figure 2. The 
curves in figure 2(a) give the maximum possible increase in critical 
stress and apply when a stiffener is used for which the stiffness is 
equal to or_ greater than the critical stiffness, that is, the value of 
stiffness beyond which a further increase does not cauee any increase 
in the buckling stress. Figure 2(b) can be used to find the increase 
in critical stress when the stiffness is below the critical value. 

In order to find the theoretical critical— stress coefficient- for 
a given curved plate- with a central chordwise stiffener (k x + Ak x ), 
first calculate the values of Z and 7 for this plate— stiff ener 
combination (where 7 is twice the ratio of stiffener stiffness to 
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plate stiffness 2 Then find the critical— stress coefficient for 

the unstiffened plate k x from figure 1. Next find the maximum possible 
increase of critical— stress coefficient Ak x due to the presence of 

a stiffener from figure 2(a) and the critical stiffness corresponding 
to this maxi TrruTn possible increase from the curve in figure 2(b). If 
the given stiffener has a stiffness equal to or greater than the critical 
stiffness , add together the stress coefficient k x for the unstiffened 
plate and the increase £k x just found from figure 2(a); the sum is the 

stress coefficient for the plate— stiffener combination. If the given 
stiffener has a stiffness less than the critical stiffness, add together 
the stress coefficient k x for the unstiffened plate and the increase £k x 

found from figure 2(b). 

The curves of figure 2 apply to plate— stiffener combinations buckling 
in one wave in either direction. As the curvature of a panel increases, 
the axial wave length tends to diminish. In order to show why the use 
of a central stiffener does not increase appreciably the critical com- 
pressive stress of the unstiffened plate when the values of Z are 
greater than those given in the preceding table, the curves of figure 3 
have been drawn for (3 = 1 to include consideration of cases of more 
than one ha lf wave in the axial direction. These curves show that for 
two half waves (node at the stiffener, li < z < 23 ) no increase in 
stress due to a stiffener occurs and that for three half waves (23 < Z < 38 ), 
only a slight increase in stress ;Ls possible and even this slight increase 
requires a stiffener of high stiffness. Only if the unstiffened panel 
buckles into one half wave in the axial direction, does it appear worth- 
while to use a single chordwise stiffener. Computations carried out for 
curved plates with the circumferential dimension 1.5 and 2 times the 
axial dimension showed the same results. In al 1 of these cases the 
minimum load was found when there was one wave in the circumferential 
direction. 

Estimation of design critical stresses .— Panels having moderate 
and large values of Z "buckle at compressive stresses considerably 
lower than the stresses predicted on the basis of the small,— deflection 
theory. In the absence of experimental data on the compressive 
buckling stress of panels with a stiffener having a flexural stiffness 
lesB than the critical stiffness, the fraction of the total possible 
increase in experimental critical stress, actually achieved by use of 
the stiffener, may be assumed to be equal to the fraction predicted by 
the theoretical solution (obtainable from figs. 1 and 2X This 
assumption permits an estimate to be made for the compressive stress 
of such a panel when the stresses for the limiting cases of an unstiffened 
panel and for a panel with a stiffener having a flexural stiffness 
equal to or greater than the critical stiffness are known. According 
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to the aforementioned assumption the critical stress cr eX p(/) for the 
stiffened panel is given by 

= CT exp^^ + ® j^exp(^cr) ~ °exp 

where 

£k x (y) 

E = ■ 

^-x^cr ) 


The increase in the theoretical critical— stress coefficients £k x 
needed to calculate E for panels having a wide range of- ratios of 
circumferential dimension to axial dimension is obtained from figure 2 . 

The stresses and °exp^ cr ^ ma 3 r regarded as the 

critical stresses of_unstiffened panels: c eX p(o) is the critical 

stress of the original panel with no stiffener and CT eX p( 7 cr ) is the 

critical stress of a panel having the radius, thickness, and circumr- 
ferential dimension but half the axial dimension of the original panel. 
Because adequate design data on the buckling stresses of curved 
rectangular panels in axial compression are not available, some method 
of approximating these stresses must-be used. Three, lower limits may 
be given for the buckling of a curved panel; they are the buckling 
stress of the corresponding flat panel which may be obtained from 
figure l(a), the buckling stress of the complete cylinder which may 
be obtained from figure 4 adapted from reference 4,_and the buckling 
stress of the long curved strip, of which the curved panel may be con- 
sidered a part, which may be obtained from figure 5 adapted from 
reference 3 . The highest -- of these three lower limits represents a 
conservative approximation to the actual buckling stress of the panel. 
By use of values for CT eX p(o), c ex p(7 cr ) 3 and E, determined in 

the manners just described, the design buckling stress cr eX p ( 7 ) may 
be determined from, equation (l). 


CONCLUSIONS 


The theoretical analysis shows that for certain curvatures and 
aspect ratios designated in this paper, a curved rectangular panel can 
be appreciably strengthened to resist additional axial compression 
without buckling by the use of a central chordwise stiffener. The 
strengthening effect decreases as the curvature increases and as the 
ratio of the circumferential dimension to the axial dimension decreases. 



BACA TN Bo. l66l 


5 


With the aid of the semi empirical results contained herein, the axial— 
compreasive buckling strength of a curved rectangular panel with a 
central chordwi.se stiffener can he estimated. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., March 11, 1948 
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APPENDIX A 
SYMBOIS 


a axial or c Ircumf erential dimension of panel, whichever 

is larger 

h axial or circumferential dimension of panel, whichever 

is smaller 


m, n,l 
Pi 1J 


integers 


r radius of curvature- of panel 

t thickness of panel 

L length of cylinder 

u displacement in axial (x— ) direction of point - in 

median surface of panel 

v displacement- in circumferential (y— ) direction of 

point““in median surface of panel 

v displacement In radial direction of point in 

median surface of panel j positive outward 

3T axial coordinate of panel 


7 

D 

E 

E 


I 

Z 


circumferential coordinate of panel 
flexural stiffness of panel per unit length 


Et~ 


,12(1 - a 2 ; 


Young* s modulus of elasticity 

Airy*s stress function for median— surface stresses 
produced "by buckle deformation ( ji> stress in 

^2 f V s r 

axial direction; — — , stress in circumferential 

\ ax 2 

direction ] 

moment of inertia of stiffener 

curvature parameter \Jl_ — p 2 f or panels or 


~ l/l — for cylinders^ 


■) 
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a kqj ®pq 


deflection coefficients in trigonometric series 


P= t 


'cr 


4 

T 


exp 


(*-t) 


critical-axial-campressive— stress coefficient 

k^D 

appearing in the formula cr-jr = — for panels 

ky-lfiV h^fc 

or cr^ = for cylinders 

l 2 t 

lateral pressure 


ratio of flexural stiffness of stiffener to half 
flexural stiffness of plate in same direction ( 2 


k 


lowest value of 7 for which a buckle node occurs 
at stiffener location 

Poisson* s ratio 

shear stress in shell 

axial— compressive stress 

circumferential-compressive stress 

experimental axial buckling stress 

operator defined in appendix C 

Dirac 5 function defined in appendix B 


Mp, - * tf*) 8 ♦ f ^ 


V - (p % 2 + i 2 ) 2 * 


z^p 8 


^ (p 2 p 2 + 1 2 ) 2 * 


- k x P 4 P 2 


v 1 *- = L_ + 2 


s? Sx 2 dy 2 dy 




inverse of defined by V - ^(y^w) =-w 
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APPENDIX B 

EQUATION OF EQUILIBRIUM OF ClfLINDRICAL SH ELL WITH STIFFENERS 


The equation of equilibrium, of a cylindrical Bhell under the action 
of shear, both direct stresses, and lateral pressure is (reference 2): 


D vK 




Sx" 


a ^ 2 


+ 2t 


a^w 
ax ay 


+ a x 


aV 


= 0 


(EL) 


In this equation each term, is a force per Unit area or pressure. The 
first term, is the restoring pressure due to the bending stiffness of 'the 
shell, the second is the restoring pressure due to the stretching 
stiffness of the shell, the third is the (negative) restoring pressure 
du.e-.to compressive stress in the axial direction, and so forth. 

The equation describing all the effects of a stiffener riveted ter - 
the shell is quite complex, since even when «.~l 1 cross-sectional deforma- 
tions are assumed zero any real stiffener has two principal bending 
stiffnesses, a torsional stiffness, a bending-torsion stiffness, and a 
stretching stiffness. In many problems it is sufficient to idealize 
the stiffener by considering it to be located along a line on the shell 
and to have only the bending stiffness that restrains radial deforma- 
tions of the shell. This Idealization is equivalent to considering the 
stiffener to be without torsional stiffness and attached to the plate 
by a frictionless bond which maintains contact but _ allows the plate to 
slide freely under the stiffener. 

If the force per unit length exerted by the stiffener on the 
panel is q and this force Is applied uniformly over the width of 
the stiffener e, the pressure becomes q/ e. A term must therefore 
be added to equation (Bl) which has the magnitude q/ e under the 
stiffener and is zero, everywhere else. As € approaches, zero this 
term, becomes simply q8 where 5 is the Dirac delta function (refer- 
ence 5)j defined by the properties 

8(t - t Q ) = 0 


wh9n 


^ 4 ^0 



NACA TN Wo. 1661 


9 


and 


bo that 




S(t - t 0 ) dt = f(t Q ) 


■when t Q Is within the internal (tq, tg) and f(t) is a continuous 
function In the neighborhood of t Q . The effect of a stiffener located 
at x = Iq would be represented by the term, q S(x — Xq) added to 

equation (Bl) and of a stiffener at y = y Q would be represented by 
the addition of this term, q S(y — y c ). 

If the stiffener runs axial ly, the beam theory gives 


n - TTT ^Nr 4. P 


(B2) 


If the stiffener runs circumferentially, the beam theory gives 


3 Nr 

<1 = El — r- 

dy 4 


+ P 


By 2 


(B3) 


where EE is the stiffness of and P the compressive force in the 
stiffener and w is the deflection of the plate at the stiffener. In 
the problem of the present paper the stiffener is located at x = ^ 

P is zero, and only axial stress is applied to the shell so that the 
equation of equilibrium becomes 


/M _ i Et „-4- a 4 w w 

\ r 2 ^ / 57 


( 


+ CT x t 


■it.l* 


■fr ■ 


* 


/;■ 

1 / 


sNr 

53 N 


= 0 


(BlO 


D v 'v 4 f ~ 


-f BJ 




u» 


l Jj 


/ 




o 
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APPENDIX C 

THEORETICAL SOLUTION POP CRITICAL AXIAL STRESS OF CfUHVED 
RECTANGULAR PANEL WITH CENTRAL CHCEDWISE STIFFENER 


Equation of equilibrium. .— The critical axial-compressive stresses 
of curved rectangular panels haYing a central chordwise stiffener of 
zero torsional and stretching stiffness placed at x = k may pe found 

"by solving the equation of equilibrium, equation (B4). Division of 
equation (Bl) by D followed by appropriate substitutions gives 


V^w + 


12Z 2 -4 a Sr 
b^ - ax 14 - 



(Cl) 


The equation of equilibrium may be represented by 

Q(w) - 0 


where Q is the operator— defined by 


Q = V k 


+ 

— F + 

b^ ax 4 


12 Z ‘ 

A 


S 8 4 

D \ 




(C2) 


(C3) 


Method of solution .— Equation (C2) may be solved by the Galerkin 
method as outlined in references 2 and 6. As suggested in reference 2 
for simply supported rectangular panels, the following series expansion 
is used for w 


v a mn sin Ip Bljl S ? 

nfc 1 n=l D 


(Ck) 


(The coordinate system used is shown in fig. 6(a).) The coefficients 
are then chosen to satisfy the equations 



NACA TIT No. l66l 


11 


Vixen the operations indicated, in equation (C5) are performed, a set of 
homogeneous linear algebraic equations in a^ is obtained -with. k x 
appearing as a parameter. The solution for the critical— axial— stress 
coefficient k^ is then found to be the mini-mum value of k x for 
■which the algebraic equations have a nonvanishing solution for a^, 
that is, for which the plate is in equilibrium in a deflected state. 

The boundary conditions implied by the method of solution are 
zero radial deflection and end moment at each edge, zero displacement 
along each edge, and free displacement normal to each edge in the 
median surface of the panel (see reference l)j that is. 


v(0,y) = 

w(b,y) 

= 0 

w(x. 

,0) 

= w(x,a) = 0 

S 2 w . . 

— - (0,y ) 


(b,y) = 0 

^w 

(*j 

,0) - (x,a) 

3x 2 

ax 2 

Sy 2 

Sy 2 

v(0,y) = • 

v(b,y) 

= 0 

u(x, 

,0) 

o 

n 

id 

•t 

*1? 

II 

& (0,y) 

_ a% 

0>,y) = o 

s 2 f 


,0) = 2_£ ( x , a ) 

dy 2 • 

dy 2 

Sx 

Sx 2 



= 0 


(06) 


Solution for circumferential dimension greater than axial dimension . 
Substitution of the expressions for Q and w given by equations (C3) 
and (Cl*-) into equation (C5) leads to the following set of algebraic 
equations : 




+ + 2| _ z£ ?VL _ _ v. .4-2 


tC (p 2 P 2 + q 2 ) 2 




_ EI A . pit 
+ 2 == sin 


Db 


a kq sin = 0 


(07) 


where 


p = 1,2, ... and q = 1,2, ... 
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which may "be written in two sets of equations as follows: 




"l>*l - ^(-D XZ .V-l) 


k-1 

2 


= 0 


(C8) 


where p = 1,3, ... and q = 1,2, ... 


and 

= .° 


(C9) 


where p = 2,b, ... and q = 1,2, ... 


For buckling across the stiffeners there must be an odd number of 
half waves in the axial direction; hence, the set of equations (C8) 
applies. Dividing through the set of equations by Np^ gives 


li P-1 oo k-1 

Vi + ^7 (-1)" 5 " 2-^ W- 1 )^ - = o 


pq 


k=l,3 


where p = 1,3, ... and q = 1,2, ... Multiplication of each equation 
by (—1) 2 and summing with respect to p gives 


P— 1 


k— 1 


2IH a pq ( " l) 2 + 7< l k 2Z 2 - a k q(-l) 2 “ 

P=l,3 Vq - p-l,3 k-1,3 


0 


(CIO) 
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■where cl = 1,2, ... Equation (CIO) reduces to 


or 


1 + ra 4 > r- - 0 

p=i, 3 %q. 


(Cll) 


7 = - 


u. 

2 


p=i>3 ^pq. 


where q. = 1,2, . . . 

It Is known that at low values of the curvature parameter Z one 
half wave (a = l) gives the minimum, "buckling load 


7 (012) 

OO 

-i- 

p=i, 3 ^Pi 

Equation (C12) gives the mini-mum value of 7 for which the plate- 
stiffener combination is in equilibrium in a deflected state at a given 
value of k x . Figures 2 and 3 present results obtained from equation (012). 

The solution to equations (C9) corresponds to no buckling across 
stiffeners and is given by 


N- 


Pl 


= 0 


where p = 2,1, ... and q. = 1,2, ... The lowest roots of these 
equations are the values of k z given in column (b) of table 1 for 
various panel shapes and stiffener stiffness. 

Solution for axial dimension greater than circumferential dimension .— 
"When the axial dimension is greater than the circumferential dimension, 
interchanging a and b in equation (Cl) is convenient in order to 
retain b as the shorter dimension 

00 00 

"N "XT „ . m . n 

-Z Z. sin — — sin 

m=l n=l a * 


w = 


(C13) 
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The coordinate system used is shown in figure 6(h). The 'use of these 
coordinates involves a slight modification of the equation of equilibrium. 

as the stiffener is now placed at x = |p. The modified equation is now 


D V*Sr + ^ + El 8 (x - () ^-r + — = 0 

r 2 & b \ V. z 


Sy* x 6x 2 


and 


Q 


.^ + 2^V+2^ + « 8 - Yx - 25_ + ^ £'S_ 

IF D V V 30* 


(014) 


(015) 


This problem is solved in a manner similar to that in the previous 
problem by a substitution of the expressions for Q and w given by 
equations (C13) and (C15) into equation (C5) (also modified by inter- 
changing a and b). The following set of algebraic equations results 




(p 2 + q 2 0 2 ) 2 


12 ZSV' , n ?_o 

+ ~V 0 “ k xP^P^ 

** ( p 2 + a 2p2 f 


E.V to T = o (0lS) 

k 

where p =1,2, ... and q = 1,2, ... 

As in the previous problem for buckling across the stiffener, an 
odd number of half waves in the axial direction is necessary and for 
no buckling across the stiffener an even number is necessary. For this 
case again at low values of Z one half wave in the circumferential 
direction gives the minimum critical stress to cause the buckling of 
the stiffener. An equation similar to equation (Cll) results, in which 
N p q is replaced by _ ' _ 


'Mp, - <p2 + q + jf Zg pV - k lP 2 p 2 


(p 2 + q 2 0 2 ) 
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and the 7 is replaced "by 7 P^. For low values of Z, q = 1 and 


7 = 


* 3 2 


P=i,3 


“Pi 


(C17) 


Equation (C17) gives the minimum value of 7 for which the plate- 
stiffener combination is in equilibrium in a deflected state at a 
given value of . k z . Few results are shown for this case as a panel 
having the axial length equal to or greater than about 1.5 times the 
circumferential length would not be significantly strengthened by 
addition of a central chordwise stiffener of zero torsional and 
stretching stiffness. (See table 1.) 

The solution which corresponds to no buckling across the stiffener 
is given by 


Mpq_ - 0 

where p = 2 ,K,6 ... and q = 1 * 2^3 ••• The lowest roots of these 
equations are the values of k z given in column (b) of table 1 . 
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IABLE 1 

CRETIGAIrTAXEAlr-CCMPEEESITB— STRESS CQESTXCXEHTS FOR BUCKLUKJ WECHODT 
AHD FOR HO BCTCKUHO ACROSS SUFFERER 


CFFERIHO ZERO TCRSIOHAIi RESTRAXRT 
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(a) Maximum possible increase. (b) Increase for given stiffness. 

Figure 2.- Increase in critical-axial -compressive-stress coefficient due to the presence 
of a central chordwise stiffener over the value for the unstiffened panel. 
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"Critical stiffness" curve 

Left of "critical stiffness" curve-, 
buckling of sheet and stiffener 

_ Right of "critical stiffness" curve - 
buckling of sheet alone 


Figure 3.- Effect of central chordwise stiffener upon critical -axial 
compressive -stress coefficients for simply supported curved 

cl 

rectangular panels of t- = 1. 
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(b) Panels with the axial dimension 
larger than the circumferential 
dimension. 


(a) Panels with the circumferential 
dimension larger than the axial 
dimension. 
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Figure 6.- Coordinate systems used in theoretical analysis. 
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